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A ROSAT HIGH RESOLUTION X-RAY IMAGE OF NGC 1068 
J. Halpern, Columbia University, New York, NY 10027 

ABSTRACT 

The soft X-ray properties of the Seyfert 2 galaxy NGC 1068 are a crucial 
test of the “hidden Seyfert 1” model. It is important to determine whether 
the soft X-rays come from the nucleus, or from a number of other possible 
regions in the circumnuclear starburst disk. We present preliminary results of a 
ROSAT HRI observation of NGC 1068 obtained during the verification phase. 

The fraction of X-rays that can be attributed to the nucleus is about 70%, so 
the “soft X-ray problem” remains. There is also significant diffuse X-ray flux on 
arcminute scales, which may be related to the “diffuse ionized medium seen in 
optical emission lines, and the highly ionized Fe Ka emission seen by BBXRT. 

INTRODUCTION 

In the hidden Seyfert 1 model of NGC 1068 (Antonucci and Miller 1985), 
the ionizing flux is seen only indirectly via electron scattering from a warm, 
ionized medium above the axis of an obscuring molecular torus. The results of 
spectropolarimetry of the Balmer lines require that the electron temperature 
in the scattering region is < 3 x 10 5 K (Miller et al. 1991, hereafter MGM). It 
is difficult, but not impossible, for soft X-rays to emerge unabsorbed through 
a gas at this low temperature (Monier and Halpern 1987). Careful fine tun- 
ing of the photoionization model is required to avoid imprinting absorption or 
emission features on the soft X-ray spectrum (MGM). This problem would be 
avoided entirely if the soft X-rays are not coming from the nucleus, a reasonable 
expectation since NGC 1068 is also a bright starburst galaxy. 

RESULTS 

A 20,000 s observation of NGC 1068 was obtained with the ROSAT HRI 
between July 24 and 26, 1990. The observed X-ray luminosity of NGC 1068 
in the 0.1 - 2.4 keV band is ~ 9 x 10 41 ergs s -1 ffor a distance of 22 Mpc), 
consistent with previous measurements. The central region of the HRI image is 
shown in Figure 1, and an X-ray contour map of the galaxy is shown in Figure 
2. There is a “nuclear” source, which is characterized by FWHM = 5'.'5 in the 
east- west direction, and 7" in the north-south direction. This is close to the 
theoretical FWHM of 5" which is expected to result when the effects of mirror 
and detector resolution, and aspect solution are combined. There is not yet 
enough calibration data available to test whether the marginal extent in the 
north-south direction is real. The fraction of X-rays from NGC 1068 that can 
be attributed to the “nucleus” is about 70%, as indicated by the X-ray flux 
within a radius of 9", which is 3 sigmas of the radial Gaussian corresponding 
to the observed FWHM. 

There is also significant diffuse X-ray flux on an arcminute scale. It covers 
much of the galaxy, and is elongated at position angle ~ 35°, which is similar 
to that of the nuclear radio structure (Wilson and Ulvestad 1983). The surface 
brightness is greatest in the northeast, which is the direction of the one-sided 
ionization cone (Pogge 1988) and the region of strong, high-excitation [Ne V] 
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FIG. 1. - The ROSAT HRI image of NGC 1068, showing the central 
16'. 7 x 16'. 7 of the field. Each dot represents one photon. Approximately 11,700 
photons are detected in a region of radius 2' centered on NGC 1068, of which 
~ 70% can be attributed to the nucleus. A stellar source is also detected 5'. 8 
south of the galaxy. 


A3426 emission (Evans and Dopita 1986; Bergeron et al. 1989). In the northeast 
direction, the diffuse X-ray emission extends as far as the large-scale 20 cm 
radio emission of the disk (Wilson and Ulvestad 1982). The diffuse luminosity 
within a 1' radius (6 kpc) is ~ 4 x 10 41 ergs s _1 when corrected for Galactic 
JVh = 3 x 10 20 cm -2 , which is comparable to that of other bright starburst 
galaxies (Ward 1988). Similar luminosities of extended X-ray emission were 
detected in Einstein HRI images of several Seyfert 1 galaxies (Elvis et al. 1990). 
In NGC 1068, about 20% of the total X-ray flux comes from a region which is 
larger than the 30" diameter starburst disk. 
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X Sky Pixel 

FIG. 2. - Contour map of the X-ray image of NGC 1068. The field is 
4' x 4'. The image has been binned into 4" x 4" pixels. The lowest contour level 
corresponds to an intensity of 4 counts per pixel, or a flux of ~ 1.5 X 10 14 ergs 
cm“ 2 s _1 . The contour levels are separated by factors of 2 in intensity. 


INTERPRETATION 

The detection of a nuclear soft X-ray source, in combination with the lack 
of oxygen features in the BBXRT spectrum (Marshall 1992), presents a signifi- 
cant challenge to the electron scattering model. The emergent X-ray spectrum 
should have either a K-edge of 0+ 7 at 870 eV (Monier and Halpem 1987), or 
a Lya recombination line at 654 eV, neither of which are seen (Marshall 1992). 
One can wriggle out of the problem by noting that the M nuclear HRI source, 
with its instrumental resolution of 500 pc, may still be outside the scattering 
region, for which the theory of MGM favors only a lower limit of 30 pc. 
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The diffuse X-ray flux is similar in extent to the “diffuse ionized medium” 
seen in the [N III and Ha imaging Fabry-Perot data of Bland-Hawthorn et al. 
(1991). Unlike H II regions, the ratio [N II] A6584/Ha is greater than 1. The 
large velocity widths of the lines in comparison with the adjacent H II regions, 
led these authors to predict that the scale height of the emitting filaments might 
be high, ~ 400 pc, and that they might be confined by, or condense out of, a 
hot phase which would be seen by ROSAT. The HRI image is substantially 
in agreement with this prediction. Since some of the X-ray flux comes from 
a region which is larger than the 30” diameter starburst disk, it is likely to 
originate in a diffuse hot medium, rather than numerous discrete supernova 
remnants or binary X-ray sources. To explain the diffuse X-ray emission as 
bremsstrahlung at T ~ 10 7 K, we require emission measure n 2 V ~ 7 x 10 64 
cm -3 , or electron density n e ~ 0.15 over the inner 6 kpc radius and 400 pc 
scale height. This medium would also be responsible for the He and H-like iron 
fines seen the in BBXRT spectrum, which requires n 2 F ~ 1.5 x 10 64 cm -3 at 
T ~ 10 7 K (Marshall 1992). It could also exist in pressure balance with optical 
filaments of density ~ 100 cm -3 and filling factor ~ 10“ 6 . 

The cooling rate of the optical filaments is estimated to be ~ 4 x 10 42 ergs 
s _1 (Sokolowski et al. 1991), while the diffuse X-ray luminosity is ten times 
smaller. Photoionization by the active nucleus is an ample source of energy for 
both the optical filaments and the X-ray emitting gas, although the ionization 
parameter may not be high enough to heat the low-density medium to 10 7 K. In 
the context of the the obscuring torus/electron scattering model, the intrinsic 
ionizing luminosity L of the nucleus is related to the observed luminosity L 0 
via L = (l/r es )(47r/fl)L 0 , where SI is the unobscured solid angle. Estimates 
of L range between 7 x 10 43 and 1 x 10 45 ergs s -1 (Monier and Halpem 1987; 
Sokolowski et al. 1991), or even higher if there is intrinsic beaming (MGM). So 
the diffuse optical filaments need only absorb ~ 1% of the nuclear luminosity to 
power the emission lines. However, the ionization parameter (£ = Tx/^r 2 ) must 
be 1000 in the X-ray emitting medium for He and H-fike iron to predominate. 
This can only occur if the nuclear luminosity is > 10 45 ergs s 1 , as estimated by 
MGM. Alternatively, shock heating by supernova remnants or energetic stellar 
winds may contribute to the energetics of the diffuse X-ray emitting gas. 
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The nature of the y-ray source ‘Geminga’ (2CG195 + 04) is a 
problem of considerable importance in high-energy astrophysics. 
First discovered in 1972 by the SAS-2 satellite 1 , Geminga emits 
virtually all its power at energies above 50 MeV, and at energies 
above 100 MeV is the second brightest source in the y-ray sky 
survey made by the Cos-B satellite 1 . It eluded identification at all 
other wavelengths until the Einstein Observatory found an unusual 
soft X-ray source, 1E0630+ 178, in its error box 3 . This source also 
has a claimed twenty-fifth magnitude optical counterpart 4-6 . This 
distinctive set of properties is reminiscent of the Vela pulsar, except 
for the absence of radio emission 7 or a synchroton nebula 3 . We 
have made a more sensitive soft X-ray observation of the Geminga 
field using Rosat, and have detected coherent pulsations from 
1E0630 + 178 at a period of 0.237 s. This result confirms 
suggestions 3-6,8 * 9 that Geminga is, like Vela, a y-ray pulsar. We 
speculate that Geminga is somewhat the older of the two. With 
this discovery we consider the mystery of Geminga largely solved. 

For Geminga, the soft X-ray band is the most promising one 
for a pulsar search because the count rate is high and the 
background is negligible. The Geminga field was observed with 
the Position Sensitive Proportional Counter (PS PC) aboard 
Rosat 10 during 1991 March 14-17. A total of 14,390 s of exposure 
were obtained during 10 satellite orbits as listed in Table 1. The 
mean count rate of 1 E0630 + 178 was 0.53 s" 1 in the 0.07-2.4-keV 
band. A total of 7,630 counts were detected within a circle of 
radius 1.5' about the source centroid, of which only ~1% are 
background. After transforming the photon arrival times to the 
barycentre of the Solar System using the accurate (~3") Einstein 
High Resolution Imager position, the search for periodicity was 
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TABLE 1 

Rosat observations of Geminga 


Date (ur) 

Start time (s) 

Duration (s) 

Counts 

1991 March 14 

76.564 

1,193 

612 

1991 March 15 

1.191 

1.437 

749 

1991 March 15 

29,362 

1,360 

703 

1991 March 15 

64,656 

1.567 

866 

1991 March 15 

76,468 

1,272 

641 

1991 March 16 

64,514 

1,690 

887 

1991 March 16 

69,999 

1,959 

1,038 

1991 March 16 

76,448 

1,272 

724 

1991 March 16 

81,745 

1,509 

786 

1991 March 17 

46,703 

1.131 

624 


begun by performing a fast Fourier transform (FFT) on the 
section of data consisting of four consecutive orbits on March 
16. The FFT covered 262,144 independent trial periods down 
to a minimum period of 0.072 s. The only significant peak in 
the power spectrum occurred at a period of 0.237097 s, for which 
the power is 17.5 times the mean. The probability that a peak 
of this strength or greater will occur by chance somewhere in 
the power spectrum is 0.007. The same section of data was then 
folded about a range of trial periods. The x 2 values of the 
resulting five-bin folded light curves were calculated with respect 
to a constant mean. The candidate period appeared with x ~ 
39.5 for 4 degrees of freedom. The probability that x greater 
than or equal to this value will occur by chance in a single trial 
is 6x 10" 8 . When multiplied by the number of trial periods, the 
chance probability is 0.016. The highest x 2 in the periodogram 
of the folded data is 54.7 (Fig. 16). The single-trial chance 
probability is less than 10~ 10 , but the effective number of periods 
searched in the epoch folding is about a factor of 5 higher than 
in the FFT. We conclude that the probability that the period is 
spurious is ^1 x 10 -4 . 

The reality of the signal at 0.237 s is further supported by its 
presence in the remaining half of the data which was not 
included in the original FFT. Figure la shows the periodogram 
for the five orbits of data obtained on March 14-15. The period 
is clearly indicated in this independent set of data, although the 
less-than-optimal sampling introduces many more aliases at the 
96-min satellite orbital period and its multiples. Figure l c is the 
result of the coherent folding of data from all 10 orbits. The 
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FIG. 1 Periodograms of the five-bin folded light 
curves. The multiple peaks are aliases of the 
96-min satellite orbital period. The best-fit period 
is 0.2370974 s±0.1 p.s. 
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highest peaks in all three panels are at the same period to within 
the errors. The symmetry of the peak structure indicates that 
the highest peak at 0.2370974 s ±0.1 pis is most likely to be the 
true period. In addition to these analyses, the period is detected 
weakly in those individual orbits having more than 1,500 s of 
exposure. There is no evidence for orbital motion or a secular 
change in period P during the span of these observations. On 
a timescale of two days, the upper limit on P is 2 x 10 -12 s s~\ 
and Au is less than 0.5kms _1 . We conclude that the 0.237-s 
period is not an instrumental artefact, because similar analyses 
on background counts collected at different positions in the 
detector do not reveal a signal. We do not detect a significant 
signal at any other period, including previously reported 
periods 11-13 near 60s, even though the PS PC is at risk of 
manufacturing spurious periods between 30 and 400 s because 
of the deliberate dithering of the satellite. 

Figure 2 shows the folded light curve in three different energy 
bands, and their sum. The pulse profile has a single broad peak, 
which is consistent with the absence of power at submultiples 
of the 0.237-s period. The pulsed fraction is 24% in the band 
0.07-0.2 keV and declines with energy. At 0.2-0.28 keV, the 
pulsed fraction is 19%, and at 0.28-1.50 keV it is 15%. These 
pulsed fractions are probably lower limits, as there is a known 
timing error in the processed data which shifts some of the 
photon arrival times randomly with respect to the phase of the 
0.237-s period. The extremely soft spectrum of the pulsations, 
and the broadness of the pulse profile, are suggestive of thermal 
emission from the surface of a neutron star. But it is difficult to 
interpret these data in terms of surface temperature and distance, 
because the X-ray spectrum in the Rosat bandpass can be 
modified from the black-body shape as a result of non-grey 
atmospheric opacity 14 . Strong L and K edges could be present. 
Rosat has recently detected soft X-ray pulsations from Vela 15 , 
also in the lowest energy channels. 

To study the X-ray spectrum, we extracted all the source 
counts within a 2' radius, and subtracted background from an 
annulus between radii of 2.5' and 4.2'. The total flux in the 
0.1-2.4-keV band is —1.5 x 10“ 12 erg cm* 2 s‘ ! , which is roughly 


FIG. 2 Pulse profile of all data between 0.07 and 1.5 keV folded into a 20-bin ► 
light curve. The top panel is the sum of the bottom three. 


consistent with the Einstein result. The X-ray spectrum is very 
soft, but it cannot be adequately fitted by a black-body or a 
power-law spectrum, as illustrated in Fig. 3. We have not yet 
explored the full range of plausible spectral models, but with 
the caveat concerning atmospheric opacity in mind, a composite 
model consisting of a black body of temperature T *= 
(3 -4) x 10 5 K plus a power-law component with energy index 
a in the range 0.75-1.75 yields an adequate fit to the data. In 
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FIG. 3 Spectral fits with single-component models, a, Power 
law; b, black body. 
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this model, the decrease of pulsed fraction with energy might 
be caused by the increasing contribution of an unpulsed syn- 
chrotron component. A column density in the range (0.5-3) x 
10 20 cm -2 is allowed. 

Observations indicating that Geminga might be a pulsar some- 
what older than Vela have been reviewed previously 3 6,s . The 
overall energy distributions of the two objects are similar at 
wavelengths from optical through y-rays. The ratio of L y / L x = 
1,000 is a likely indicator of a pulsar emission mechanism. The 
ratio LJ L opl ~ 1,800, and the soft X-ray spectrum, also suggest 
an isolated neutron star. The discovery of a short pulse period 
fulfils this expectation. A detailed theoretical model for y-ray 
pulsars was developed by Ruderman and Cheng 9 , who proposed 
an evolutionary scheme in which a short-period y-ray pulsar 
evolves from a Crab-like stage, through that of Vela, and finally 
reaches that of Geminga. In the regime of Vela-like pulsars, 
y-rays are produced in the outer magnetosphere as synchrotron 
radiation from electron-positron pairs. As it slows down, such 
a pulsar uses an ever increasing fraction of its spin-down power 
in the production of an e* wind and associated y-ray emission 
until almost the full spin-down power, -4 x 10 36 erg s _l , goes 
into this mechanism. As the period increases further, the currents 
necessary to sustain pair production are quenched and the y-ray 
emission ceases. They predicted that Geminga is near the turnoff 
period, which for Vela itself would be — 0.13 s, and postulated 
that ordinary pulsars avoid this fate because they have more 
highly curved magnetic field lines near the polar cap, which can 
maintain the potential drop required for pair production. 

The theory was generalized by Chen and Ruderman 16 , who 
argued that there is actually a death line for Vela-like (outer 
gap) pulsars given by 5 log B — 12 log P = 72.2, where B is the 
surface magnetic field. If Geminga is to fall near this line, it 
must have B «= 1 x 10 13 G, only a factor of 3 greater than that of 
Vela. In this case, the total spin-down power, E ~ B 2 R 6 D, A /c\ 
would be -1.8 x 10 36 ergs* 1 . As the observed y-ray flux 1 is at 
least 2.4 x 10 -9 erg cm* 2 s" 1 , an upper limit to the distance is 
2,500 ( B/ 10 13 ) pc for the luminosity not to exceed the estimated 
spin-down power. (The column density derived from the X-ray 
spectrum almost certainly restricts the distance to less than 
500 pc.) The predicted magnetic field would result in a spin- 
down rate of 4x 10" 13 (fi/10 13 ) 2 s s" 1 , significantly below the 
upper limit of 2 x 10* 12 established by this observation. A second 
observation could easily determine P, and thus the magnetic 
field. The fact that Geminga is not detected in high-energy X-rays 


or low-energy y-rays 17 is consistent with the outer gap model. 
The low-energy cutoff of the spectrum at a frequency o> c is 
determined by the energy of the e* pairs that leave the magneto- 
sphere before losing their energy to synchrotron radiation. 
Because a > c ocQ~ 2 B~\ scaling from the Vela break at M MeV 
(ref. 7) predicts a break at -25 MeV for Geminga. Other pulsars 
that may be progenitors of the Geminga-like y-ray pulsars are 
PSR1509-58 (ref. 18) with P = 0.150s and PSR1706-44 with 
P = 0.102 s. These both fall in the region of Vela-like pulsars 16 , 
and were recently detected 19,20 by the Compton Gamma Ray 
Observatory (CGRO). 

Our results thus largely confirm the expectation that Geminga 
is a pulsar similar to Vela. The absence of a synchrotron nebula 
around Geminga remains a fundamental difference, which is 
perhaps related to the fact that Geminga has the longest period 
of the y-ray pulsars. We suggest that many of the other un- 
identified high-energy y-ray sources are spinning neutron stars 
which are ‘radio quiet*, or whose radio beams do not intersect 
the Earth. Rosat should be able to identify y-ray pulsars by 
finding the periods of sources detected by CGRO in the galactic 
plane. ^ 
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